Radiative Corrections to Neutrino Deep 
Inelastic Scattering Revisited 



A.B. Arbuzov a , D.Yu. Bardin b , and L.V. Kalinovskaya b 

^Bogoliubov Laboratory of Theoretical Physics, 
JINR, Dubna, 141980 Russia 

b Dzhelepov Laboratory of Nuclear Problems, 
JINR, Dubna, 141980 Russia 



Abstract 

Radiative corrections to neutrino deep inelastic scattering are revisited. One-loop 
electroweak corrections are re-calculated within the automatic SANC system. Terms 
with mass singularities are treated including higher order leading logarithmic cor- 
rections. Scheme dependence of corrections due to weak interactions is investigated. 
The results are implemented into the data analysis of the NOMAD experiment. The 
present theoretical accuracy in description of the process is discussed. 
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1 Introduction 



Modern experiments such as NOMAD [1], NuTeV [2], and CHORUS [3] made a serious 
step forward in studies of neutrino deep inelastic scattering. Their precision measure- 
ments made it necessary to update the accuracy level of the theoretical description of the 
process. Important ingredients for an advanced precision in the theoretical predictions 
is the calculation of the relevant radiative corrections (RC). 

Our study is motivated by the request from the NOMAD experiment. Electroweak (EW) 
radiative corrections to neutrino-nucleon scattering should have been implemented into 
the general Monte Carlo system for the experimental data analysis. Certain experimen- 
tal conditions of particle registration and event selection should have been taken into 
account. In order to make the relevant subroutine describing the corrections fast, we 
had to look for an analytical answer (without numerical integrations). 
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In the present work we reproduce most of the results of unpublished communication [4] . 
Contrary to the earlier calculation, we used the modern technique of automatic calcula- 
tions within the SANC system [5,6,7,8] developed for Support of Analytic and Numeric 
calculations for experiments at Colliders. Besides the one-loop calculation, we consider 
certain contributions of higher orders and discuss the theoretical uncertainty due to 
unknown electroweak (EW) corrections for the case of the concrete experimental study. 

The paper is organized as follows. In the next section we define the notation and present 
the Born level distributions. Then we consider different sources of radiative corrections. 
Special subsections are devoted to quark and muon mass singularities. Numerical results 
and comparisons are presented in Sect. 4. Possible applications of our formulae are 
discussed in Conclusions. 



2 The Born Cross Section 



We will consider the process of deep inelastic scattering (DIS) in the framework of the 
quark-parton model, assuming that we are in the proper kinematical region (Q 2 ^> 

Here we list the Born level cross sections of neutrino-quark interaction weighted by the 
quark density functions. For the charge current (CC) scattering processes 



v(fa) + QiiPi) -> l (k 2 ) + qf(p2) and v + q { -> l + + q f (1) 
we have 
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= v°cc(i-y) 2 , (2) 



r°cc = \Vif\ 2 -?3 , % „ w ~ f^Q 2 ), (3) 



where M w is the W-boson mass; s is the center-of-mass energy of the neutrino-quark 
system squared; |V^/| is the element of the Cabibbo-Kobayashi-Maskawa quark mixing 
matrix; fi(x, Q 2 ) is the density function of the initial quark in the given nucleon. The 
kinematics 1 is described by the Bjorken variables: 



Q 2 A2 / \2 Q 2 

y = —, Q =-(p2-Pi), x = — : 

s yS 
s = (h+ Pl ) 2 ^xS, S=(h + P) 2 , (4) 



1 We use the (+, —,—,—) metrics, p = (p 



2 



where P is the initial nucleoli momentum. We will use mi^i (Qi,2,i) for the masses 
(charges) of the initial quark, the final quark, and the muon. 

For the neutral current (NC) scattering processes 



v(h) + ft(pi) -> v(h) + 5/0*0 and v + q t -> 1/ + <?/ (5) 
we have 



d 2 ^^™ 0/2,2/1 \2 2m 2 y 
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r' 2 A/f 4 

^c = ^— /^Q 2 ), (8) 



where 



9l = ~ + \Qi\ sin 2 6 w , g R = \Qi\sm 2 w , (9) 
M z is the Z-boson mass; and 9 W is the weak mixing angle. 



3 Radiative Corrections 



The radiatively corrected neutrino DIS cross section can be represented as the sum of 
the Born distribution with the contributions due to virtual loop diagrams (Virt), soft 
photon emission (Soft), and hard photon emission (Hard): 



J 2 Corr. j 2 Born 

U °' U ' • srVirt , rSoft , rHard 



sV=ijVl 1+i ' +< < +s >> < 10 > 

where the index % denotes the type of the process under consideration [yq CC, vq NC 
and so on). This formula assumes calculation of the order O (a) RC. Certain higher 
order corrections will be added below as well. We assume also that the momentum 
transfer square is small compared with the W^-boson mass: corrections of the order 
aQ 2 /M^ are omitted from the calculation. 
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3.1 Virtual Corrections 



Contributions due to virtual (one-loop) EW corrections are re-calculated by means of 
the SANC system. All the contributing one-loop diagrams were calculated in the frame 
of the Standard Model. Omitting the terms suppressed by the Q 2 /M^ ratio allowed us 
to get the short analytical answer. 

The virtual correction to the CC neutrino-quark scattering reads 



rVirt _ a 

°uCC — ~ 
71 
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Q 2 



Q 2 ln 2 ^ I + (l + g 1 ) 2 ln 2 4 + ln2 4 + ln ' • 

///.] mf 

s 



mf 
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ml 



Q,{2 + gO In 2 y - 2Q 1 ln 2 (l - y) - Q 1 (6 + 5Q0 lny 



+6Q! ln(l -y)- 4(2 + Qx) 2 C 2 + 1 + Qi + 8Qj 



|d + Q.)ln^}, 



7T 



C 2 = Li 2 (1) = — , 



Li. 



<*> = -/ 



ln(l — xy) 



where A is the auxiliary photon mass, A <C m\^i- Note that in the above expression we 
used the explicit value Qi = — 1 and eliminated the final state quark charge by applying 
the charge conservation law, Q 2 = Qi — Q\. 



For the CC antineutrino-quark scattering we get 
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For the pure QED part of the virtual corrections to the NC neutrino (and antineutrino) 
scattering we have 



'NC 
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■ In 2 ^ • - In 2 Q ~ 
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1 

4 



Q 

m l 
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777t 
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(13) 



The weak part of the virtual correction to the NC case is included in the definition of 
the effective electroweak couplings: 



9l(r) -> 9l(R), 9l = p( - 2 + IQil^sin 2 ^, = sin 2 6» w , (14) 

where p and k are the electroweak form factors, 



3 f 1 , (Ml\ M 2 H 
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+ 25 / (Q 2 ;m / ,m / )-n^(Q 2 ), 



(15) 



(16) 



where /]• , c/, Qf, Vf and m/ are weak isospin, color factor (1 for leptons, 3 for quarks), 
charge, vector coupling and mass of a fermion; notation for W, Z, H, t quark masses are 
evident; s w and c w are sine and cosine of weak mixing angle and = 1 — M^/Mj; 
Ap is Veltman parameter, 



Ap = -L [E«) - S(M 2 ) 

1V1 W 



(17) 



function Bf(Q 2 ;mf,mf) is a useful combination of the finite parts of the standard 
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Passarino-Veltman functions (without UV pole 1/e) 



B f (Q 2 ;m f ,m f ) = 2 B 2 i(Q 2 ;m f ,m f ) + B 1 ( y Q 2 ;m f ,m f ) ; 



(18) 




(19) 



/ 



Note, that the three last quantities are taken at /i = M w , with fi being the t'Hooft scale. 
More about the calculation and renormalization scheme applied here can be found in 
the book [9]. 

3.2 Soft Photon Radiation 

Emission of a soft photon in neutrino DIS can be described in the standard way by the 
accompanying radiation factors: 



We consider the problem in the rest reference frame of the final quark, p 2 = 0, which is 
equivalent (in the soft photon limit) to the so-called /^-reference frame, used below in 
the evaluation of the hard photon contribution. The soft photon energy, cu, is limited 
by the parameter ou, which is assumed to be small compared to the large energy scale: 
oj <C \FQ 1 ■ List of the relevant integrals is given in Appendix A. For the charged current 
we get 
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Soft corrections for the NC case are given by 



rSoft 
°NC 
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(24) 



The above formula is valid both for the neutrino-quark and antineutrino-quark NC 
scattering processes. 



3.3 Hard Photon Radiation 



In the case of hard photon emission, like in the process 



we 



f(fci) + ft(Pi) Z (fe) + 9/(p2) + 7W> 
have to extend the list of kinematical variables: 



(25) 
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p 2 =p 2 + k, Q 2 = -2 Pl (p 2 + k), M 2 = -{p 2 + k) 2 . 



(26) 



Note, that in the case without real photon radiation Q 2 = Q 2 . Therefore in what follows 
we can use Q 2 instead of Q 2 . The calculations of the hard photon contribution were 
evaluated also in the environment of the SANC system. 

The angular phase space of the hard photon is dflR = d cos 0r<1(Pr in the i?-system of 
reference, where k + p 2 = 0. 

For the charged current uq scattering we have 
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Hard photon contribution to the uq CC scattering reads 
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- 7 -(l - yfQ\) lny + Q - \y + y 2 - |(1 - y)^) ln(l - y) 
_(1 _ y )2 (i _ 2 Q 1 ) Li 2 (y) - (1 - y) 2 (l - - \ + yQi 
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For the neutral current scattering we have 
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//t 

17. 2 1 2 5 15 1 2 7 
-12 ln ^-2 ln ^ + 4 ln ^i + 2 ln ^"I 1 ^ 
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The corresponding correction to the NC antineutrino-quark scattering, 5?^, can be 
obtained from the above equation by the substitution g\ <-> y^. 



3.4 Quark Mass Singularities 



In the sum of the soft and hard photonic corrections, the auxiliary parameter lD (soft- 
hard separator) cancels out. The infrared singular terms (with logarithm of A) cancel 
out in the sum of the virtual and soft contributions. Moreover in the total sum, the 
large logarithms (mass singularities) with the final state quark mass m 2 disappear in 
accordance with the Kinoshita-Lee-Nauenberg theorem [10,11]. 

Large logarithms containing the initial quark mass, \n(Q 2 /m\), remain in the sum of 
all contributions. But these logs have been already effectively taken into account in 
the parton density functions (PDFs). In fact, QED radiative corrections to the quark 
line are usually not taken into account in procedures of PDF extraction. Moreover, the 
leading log behaviors of the QED and QCD DGLAP evolution of PDFs in the leading 
order are proportional to each other. So one gets evolution of PDFs with an effective 
coupling constant 



a 



eff 



a s + 



Ql 

a 



-a, 



(30) 



where a s is the strong coupling constant, and Cp is the QCD color factor. The nontrivial 
difference between the QED evolution and the QCD one starts to appear from higher 
orders, and the corresponding numerical effect is small compared to the remaining QCD 
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uncertainties in PDFs [12,13,14]. See also Ref. [15] for a discussion with respect to the 
process under consideration. 

The best approach to the whole problem would be to re-analyze all the experimental 
DIS data taking into account QED corrections to the quark line at least at the next- 
to-leading order. But for the present task we can limit ourselves with application of the 
MS subtraction scheme [16] to the QED part of the radiative corrections for the process 
under consideration. This leads to a shift of the initial quark mass singularities (with a 
certain constant term) from our result into the corresponding quark density function. 
The latter should be also taken in the MS scheme 2 . In fact, using the initial condition 
for the non-singlet NLO QED quark structure function (which coincides with the QCD 
one with the trivial substitution Cpa s — > Q 2 a, see Ref. [17]), we get the following 
expression for the terms to be subtracted from the full calculation with massive quarks: 



Qi 



a 
27 
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dx x 



1 + x 2 
1-x 



In 
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mf 



- 1 - 21n(l -x) 



+ 



^ttV 3 m? 9 J 



(31) 



3.5 The Muon Mass Singularity 



In the case of CC scattering, the large logarithms singular in the limit mi — > re- 
main in the final answer. These terms are in agreement with the prediction of the 
renormalization group approach. In fact, they can be described by the electron (muon) 
fragmentation (structure) function approach [18,19,20]: 



da" = X)d*i® Q 2 ), (32) 

j 

where d&j is the differential cross section of neutrino DIS with production of particle 
j {j = Hi 1) etc.); D(£,Q 2 ) describes the probability to find a muon with the relative 
energy fraction £ in particle j; sign <g> stands for the convolution operation. 

For our purposes it is enough to consider only the leading log approximation in the 
O (cxL) and O (a 2 L 2 ) orders including the contribution of electron-positron pairs. Under 
these conditions we can write the QED fragmentation function in the form: 



2 A similar consideration can be performed in the DIS and other schemes as well. 
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(33) 



where L is the so-called large logarithm, and -P^(£) is the lowest order non-singlet 
splitting function: 



P (o) (0 = 



i+e 



(34) 



For the definition of the plus prescription and the convolution operation see for instance 
Refs. [21,20]. 

Differential cross section d&j (for j = /x) is the Born level distributions (2) with a shifted 
value of variable y: 



y 



(35) 



which provides the proper value of variable y after the fragmentation stage under the 
adopted condition for particle registration. Convolution with the splitting function 
gives the following first order leading logarithmic corrections: 



&c L = ^L(2]ny-Hl-y) + ^-y), 



uCC (l-y) 2 2tt V v ; 1 



3 

y 2 



(36) 
(37) 



We will see that the above expressions reproduce the the leading log part of the complete 
result for the O (a) correction in Eqs. (41,42). 

By convolution with the O (a 2 ) term from Eq. (33), we get the second order leading 
log corrections to the CC neutrino DIS: 



€i L c L 



a 



2 V2tt 



L 



4 In 2 y - 4 In y ln(l - y) + X - In 2 ( 1 - y) - 4( 2 



+(6 -Ay) lny+(3y-4) ln(l - y) + 
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1 a 
3 27" 
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>vCC 



1 



1 / 2 



[i-y? 2 V2tt 



1 



4(1 - yf ( Li 2 (1 - y) + In 2 y + - ln 2 (l - y) 



-■\> ) (6 1.7) hi // - - ( y - -) ln(l - y) + - - 2y 



1 a 
3 



, r;(l)" 
+ - —LOdcc ■ 



(38) 



(39) 



Argument of the large logarithm depends on the energy scale of the process under con- 
sideration and the muon mass. We have two scales: s and Q 2 . In the actual calculations 
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the logarithm with the muon mass singularity arises in the form L = ln(s/m 2 ). So we 
choose s, while the difference between the two possibilities appears in our case only in 
terms of the order O (a 2 L), which are omitted now in any case. 

Using the formalism of the fragmentation function approach [22] we can get the next- 
to-leading corrections 3 of the order O (a 2 L). But, as will be seen from our numerical 
estimates below, they are small compared with the requested precision tag. 



4 Numerical Results and Conclusions 



Summing up all the different RC contributions considered above, and applying the 
MS subtraction of the initial state quark mass singularity we arrive at the result for the 
corrections to neutrino-quark cross sections: 



j2_Corr. j2_Born 

d "' ~ ia - 'l + * j + 4 2)11 -<W), (40) 



dx dy dx dy 
where 5^ LL vanishes in the NC case. 
For neutrino-quark CC scattering we have 

*»cc = I { ip z + (i - \y - ^ Mi - v) + in y) in A> 

1 3 11 

+ 2 ln ( 1 - y) ln y - 2 Li2 ^ ~ 2 ln2 y ~ 2 ln2(yl ~ y>> 

7 1 5 1 

+ (1 - y) ln(l - y) - - lny + -y lny + -y + - + 2( 2 

+ Ql ( 3 ~l ln Jp +(2 ' ln(1 ~ V) lny ~ 2Li2 
2, Q 2 23 5 1 2 .\1 

where we used the charge conservation law Qi = Q2 + Q1 and the explicit value Qi = — 1. 
For antineutrino-quark CC scattering we have 

*"» = JTW -«{ ( (1 " ^ ^ " I + i) ln 4 + (1 " s)2 ( 2C2 " Ll2 fe) 

+ ln(l - y) lny - 1 In 2 y - ln 2 (l - y)) + (-^ + ^y - y 2 ) ln 



See Ref. [23] for an analogous calculation. 
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4 + v + Qi 



1 5 7 2 3 2 s(l-y) 

2 "a* +4* -^-V) 



+ (l-y) 2 (2Li 2 (j/)+lii(l-y)lii J /-C 2 

-(i-y)C2- Ig y -gy + ^ 



+ Q 2 



2 n v 2l Q 2 



(42) 



For the case of NC neutrino scattering we get 



a 



-Qi 



Q 



l 



4-6(i-,)') +S S(>S + s( 1 -'>' + 5( 1 -') 



^-C2)]+»!+s|(l-B) 2 }-l. 



(43) 



The expression for the one-loop radiative correction to the antineutrino NC process 
can be received from the above one using the substitution: 



<Wc = $uNc(9l <-> 9r, 9l <-> 9r)- 



(44) 



Formulae (41,42,43) completely agree with the ones derived in Ref. [4]. Moreover, our 
formulae for the CC scattering case agree with the calculations presented in Ref. [24], 
where the explicit expressions were given for vd and vu CC scattering channels. For 
the comparison, one should subtract from their formulae for g"(y, S) and g"(y, S) the 
quantities C and (1 — y) 2 C, respectively, 



C 



M 3 
ln w 



1 

/iMJ + 4' 



(45) 



The subtraction reflects the choice of the MS renormalization scheme. 

Let us consider numerical results obtained for the following conditions: the fixed neu- 
trino energy E v = 80 GeV; isoscalar nuclear target; cut on the energy of the final state 
hadronic system -Ehadr > 10 GeV. We used the CTEQ4L set [25] of parton density func- 
tions. In Table 1 we present numerical values and the absolute shifts due to radiative 
corrections of the quantities [26], constructed from the cross sections of neutrino DIS, 



" a» cc (^N - fj,-X) 



(46) 
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(ISCH.ILLA), 

EW scheme 


R u 


5R U NC 


SR cc 


A^sin 2 6 W 


R- 


A" sin 2 9 W 


(0,1), G F 


0.31006 


-0.00291 


-0.02147 


-0.01130 


0.27094 


-0.00737 


(1,1), G F 


0.31063 


0.00071 


-0.02327 


-0.01044 


0.27195 


-0.00636 


(1,2), G F 


0.31067 


0.00071 


-0.02315 


-0.01039 


0.27196 


-0.00637 


(1,2), a(0) 


0.31080 


-0.05816 


0.03743 


-0.01020 


0.27209 


-0.00623 



Table 1 

Effect of RC on R v , R~ , and sin 9 W values in different approximations. 



R- 



H-X) - a^ci^N 



_Corr. Born 
uNC ~ a uNC 
_Born 
vNC 



,-Corr 



Born 



SRcc — 



°vCC ~ °v CC 
_Born 
°vCC 



(47) 
(48) 



The Born level values are R V Q = 0.31764 and Rq = 0.27831. In the computations we 
used the following values of constants and parameters taken from [27] (the same as the 
ones given by Eq. (4.6) of Ref. [15]). Flag ISCH defines the choice of the initial quark 
mass singularity treatment: ISCH=0 means no any subtraction (as has been done in 
Ref. [4]), and ISCH=1 corresponds to the MS subtraction scheme as discussed above. 
The order of the series in the leading large logarithms with the lepton mass singularity 
is governed by flag ILLA: O(aL) for ILLA=1 and 0(a 2 L 2 ) for ILLA=2. 

For the illustrations we used the simple tree-level relations between the shifts of R v ~ 
and sin 2 6 W : 



5R V = (l - — sin 2 d w J& v sin 2 9 W , 5R~ = A" sin 2 9 W , (49) 

where 5R U and 5R~ are the differences between the corrected quantities and the tree- 
level ones. In a real case of experimental data analysis the contribution of radiative 
corrections should be estimated within a general fit including all other effects like de- 
tector efficiencies, nuclear shadowing etc. 

One can see by comparing the first two lines in Table 1, that the MS subtraction of the 
initial quark mass singularity gives a numerically important contribution. We considered 
two electroweak schemes: the Gp erm i and a(0) (see Ref. [15] and references therein for 
the definition). We found that the EW scheme dependence is visible, but less than the 
one observed in Ref. [15]. 

One can compare our numbers with the corresponding results given in Table 1 of 
Ref. [15] for the case of the energy cut E^^ phot > 10 GeV. The difference between the 
two calculations appeared to be not small. Partially this is related to different treat- 
ments of the effect of the scaling violation in the parton density functions. A more 
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detailed comparison is in progress and will be published elsewhere. 

In Figures 1,2 we present the relative value of the sum of radiative corrections, 

for two channels (neutrino— up quark NC scattering and neutrino— down quark CC 
scattering) as a function of y for three fixed values of x. The neutrino beam energy and 
other parameters are the same as the ones defined for Table 1. One can see that the 
behavior of the corrections is rather smooth. Note that in the case of a different choice 
of variables, when Q 2 is defined from the observation of the outgoing muon, the size of 
the corrections for the CC case becomes several times larger 4 . 

-1.5 



-2 

o 

| -2.5 

I/O 

-3 



-3.5 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

y 

Fig. 1. Relative effect of radiative corrections to v — u NC scattering as a function of y for 
three fixed values of x. 
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y 

Fig. 2. Relative effect of radiative corrections to v — d CC scattering as a function of y for 
three fixed values of x. 

Let us consider the sources of the theoretical uncertainty, related to the incomplete 
knowledge of the radiative corrections. The sources are listed in Table 2. First, we have 
observed already the electroweak scheme dependence. From the general point of view, 
the GFermi scheme looks preferable for our problem. By trying the other EW scheme, 
we get an estimate of how large can be unknown higher order electroweak corrections. 

4 The same effect is well known in the conventional charged lepton deep inelastic scattering. 
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Source of uncertainty 


Estimated value 




for R" 


for A u sin 2 8 W 


for R- 


higher order EW RC 
higher order QED RC 
0(aa s ) 


1.5 • KT 4 
0.5 • 10~ 4 
2.2 • 10~ 4 


2.2 • 10~ 4 
0.7- 10~ 4 

3.3 • 10~ 4 


1.2 • 10~ 4 
0.2 • 10~ 4 
1.5- 10~ 4 


all combined 


2.7- 10~ 4 


4.0 • 10" 4 


2.0- 10~ 4 



Table 2 



Estimates of different contributions to the theoretical uncertainty. 

The pure QED corrections 5 are in our case not large, and by looking at the leading 
second order logarithmic contribution, we can put a limit on unknown higher order 
QED contributions. As the main source of the theoretical uncertainty we consider the 
contribution of the radiative corrections of the order O (aa s ). They can appear as in the 
loop insertions into the PU-boson propagator as well as in amplitudes where photonic 
and gluonic lines appear simultaneously. In our approach, we separated the EW RC 
from the QCD ones. After all, we receive a direct product of QCD effects (see i.e. 
Ref. [28]) and EW corrections, which doesn't cover the full O (aa s ) answer. Moreover, 
as can be seen from our calculations, we always consider Q 2 to be the argument of the 
partonic density function. But in the case of hard photon radiation off lepton, Q 2 doesn't 
coincide with the square of the "true" hadronic momentum transfer. This leads to a 
certain effect in O (aa s ). To estimate the corresponding uncertainty we varied the value 
of Q 2 in the argument of the PDFs. As the result of our estimates we found that our 
result for radiative corrections can receive up to ±3% of a relative shift. Having a unique 
value of uncertainty for the functions of two variables seems to be reasonable, because 
the functions are relatively flat, and because, in the derivation of the number we looked 
both at the variations at the differential level and at the ones for the integrated cross 
sections. Note that we estimated the uncertainties in the form of maximal variations. 
A work is going on to extend our consideration to the case of other variables and to 
avoid the effect of improper argument in PDFs, mentioned above. 

It is worth to note, that radiative corrections to the same series of processes calculated 
in a different set of variables can have a completely different behavior. For example, 
in the variables adopted by the NuTeV experiment, the size of corrections to the CC 
scattering (where the energy and angle of the outgoing charge lepton can be measured) 
are much bigger than the corrections presented here. In the same manner, estimates 
of the theoretical uncertainties should be considered taking into account the choice of 
variables and, may be, other relevant experimental conditions. Nevertheless we agree 
with the authors of Ref. [15], who claimed that the higher order contributions to the 
theoretical uncertainty seem to have been underestimated by the NuTeV experiment [2] . 

Our analytical results are combined into a FORTRAN code, which can be directly used 



5 Unambiguous separation of EW and QED corrections in the CC channel is not possible. 
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for experimental data analysis. As can be seen from the numerical estimates the effect 
of radiative corrections is greater than the present experimental uncertainty, see e.g. 
the NuTeV [2] result: 



sin 2 9 W = 0.2277 ± 0.0013(stat.) ± 0. 0009 (syst.). 



(50) 



That makes it important to take RC into account in the proper way. Estimates of the 
theoretical uncertainty becomes also relevant for the derivation of the resulting error 
in the analysis of experimental data in the present and future precision experiments on 
neutrino DIS. 
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Appendix A 

List of Integrals for Soft Photon RC 

All relevant integrals are defined in the reference frame p 2 = 0. They are: 




(A.l) 



(A.2) 



(A.3) 



/ 
/ 
/ 



d 3 k 2pip 2 
io {p\k)(p 2 k) 

d 3 k 2p 2 k 2 
uj (p 2 k)(k 2 k) 

d 3 k 2pik 2 
to (pik)(k 2 k) 




2n (21n^L 1 -iL 2 -2C(2)) ) 
2n (21n^-iL 2 -2C(2)), 



(A.4) 



(A.5) 




(A.6) 
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2 1 2 ' 1 2 2 2 ' 

mf mf wij mf mfm^ 

w = 2pik 2 w s - Q 2 , Q 2 = 2m 2 p?, 

1 + c l + cos(piA; 2 ) , m|(s - Q 2 ) 

= = 1 - pa 1. 

2 2 sQ2 

Terms which are suppressed by the ratio of a quark or muon mass to a large energy 
scale (beam energy or momentum transferred) are neglected. 
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